Role of the mean curvature in the geometry of 
magnetic confinement configurations 

A. A. Skovoroda * LA. Taimanov t 



Abstract 

Examples are presented of how the geometric notion of the mean 
curvature is used for general magnetic field configurations and mag- 
netic surfaces. It is shown that the mean magnetic curvature is related 
to the variation of the absolute value of the magnetic field along its 
lines. Magnetic surfaces of constant mean curvature are optimum 
for plasma confinement in multimirror open confinement systems and 
rippled tori. 

Introduction 

The mean curvature is one of the basic notions in the geometry of 
surfaces and vector fields [1, 2], while the geometry of magnetic fields, 
as applied to the problems of fusion plasma confinement, does not 
in fact invoke this notion (as is confirmed in [3]). The object of the 
present methodological note is to bridge this gap. We begin by giv- 
ing the definitions and main results on the differential geometry of 
surfaces in three-dimensional Euclidean space and of vector fields and 
then present various applications of the notion of mean curvature in 
magnetic confinement systems. 
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1 Mean curvature of surfaces and vec- 



tor fields 



The mean curvature of a surface, Hs, is a local quantity denned by 
the formula 



where fc m i n and & max are the minimum and maximum curvatures of 
the lines of intersection of the surface with mutually perpendicular 
planes passing through the normal to the surface at its given point. 
For e - (V x e) = 0, draw a family of surface orthogonal to a unit vector 
field e (such that |e| = 1). The mean curvature of these surfaces is 
called the mean curvature of the vector field e. In the general case 
e • (V x e) / and the mean curvature H of a vector field e is defined 
by the formula [2] 



For a unit vector e = n along the normal n to a surface, the general 
definition of the curvature H coincides with that of Hs (see the def- 
inition (1)). The surface of zero mean curvature, Hs = 0, is called a 
minimal surface. Examples of minimal surfaces are given by surfaces 
of minimum area S with a fixed boundary (soap films). Closed mini- 
mal surfaces do not exist. Among surfaces of constant mean curvature 
Hs = const are soap films between media at different pressures p, in 
which case the mean curvature is just the pressure difference. An- 
other example of the surfaces of constant mean curvature Hs = const 
are those among the surfaces bounding regions of given volume V that 
have a minimum area S. Such surfaces are called isoperimetric profiles 
and are perfect spheres of constant curvature. The Aleksandrov theo- 
rem [4] states that such spheres are the only embedded (or nested, i.e. 
non-self-intersecting) closed surfaces of constant mean curvature. A 
consequence of the theorem is, in particular, that embedded (nested) 
tori of constant mean curvature do not exist. 

The statement of this general theorem can be simply verified for 
surfaces of revolution. Assuming that the axis of revolution is the 
Z axis of a cylindrical coordinate system, we specify the shape of a 
surface of revolution by the equation <E> = z — z{r). Substituting the 
normal to this surface, n = V<J>/| V3>|, into (2), we obtain the equation 
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where the prime denotes the derivative with respect to the radius r. 
For a constant mean curvature, Eq. (3) is integrable: 

Hsr2+0 =vms' (4) 

where C is a constant of integration. Equations (3) and (4) imply that 
a closed plane curve having two or more points where z' = at r > 
does not exist. Consequently, a torus of constant mean curvature is 
impossible. Unfortunately, in considering an axisymmetric example in 
[5], sad mistakes were made that led to an erroneous conclusion about 
the existence a torus of constant mean curvature. 

Equation (4) shows that there may be two points where z' = oo 
at r > 0. The generating contours of these surfaces of revolution are 
described by a focus of a hyperbola 1 or an ellipse (see Fig. 1) rolled 
along the straight axis of revolution [6]. 




Figure 1: Rolling of hyperbola and ellipse. 



A cylinder is a limiting case of rolling of a circle. This demonstrates 
the existence of surfaces with H$ = const that are periodic along the 
axis of revolution (see Fig. 2). 

In the theory of surfaces, the Willmore functional determined by 
an integral of the square of the mean curvature over a surface plays 
an important role. For closed surfaces, this functional is conformally 
invariant: in a conformal mapping of a three-dimensional space onto 
itself, the values of the Willmore functional on a surface and on its 
transform coincide [7]. This functional not only plays a key role in 
the Weierstrass representation of surfaces [7] but has recently come 
into use in biophysics and colloidal chemistry — disciplines in which it 
is known as the Helfrich functional [8]. The critical points of the Hel- 
frich functional generate Willmore surfaces [7] . In contrast to minimal 

A rolling of a hyperbola determines one period of continuous curve in Fig. 1. 
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surfaces, there are many examples of closed Willmore surfaces, includ- 
ing embedded (non-self-intersecting) ones. Among them are, e.g., all 
spheres of constant curvature and also a Clifford torus — a surface of 
revolution generated by revolving such circle about its axis that the 
ratio of the distance R between the center of the circle and the axis 
of revolution to the circle radius r is y/2. 

2 Mean curvature of a magnetic field 

According to (2), a magnetic field B = B ■ b, where |b| = 1 and 
B = |B| is the absolute value of fiel, can be characterized by the mean 
curvature 



The variation of the absolute value of the magnetic field along its lines 
is an important parameter of plasma magnetic confinement systems 
and is often used in the geometry of magnetic fields [3, 9]. For a 
vacuum magnetic field B = — V99, where <p is the scalar magnetic po- 
tential, the mean curvature coincides with that of an equipotential 
surface (p = const for which the vector b is a unit normal vector. Note 
that, by (2), the expression (5) for the mean curvature is also valid 
for VxB/0. 

At the extreme points of the absolute value of the magnetic field 
along its lines, we have Hf, = 0. In the so-called isodynamic toroidal 
configurations revealed by D. Palumbo [10], the absolute value of the 
magnetic field is constant on the nested magnetic surfaces, B = const, 
so the mean curvature is zero: Hf, = 0, over the entire confinement 
region. Such isodynamic toroidal configurations are possible only in 



Figure 2: Nodoid and unduloid. 




(5) 
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the presence of the discharge current [11]. Without the discharge 
current, i.e., in vacuum, the equipotential surfaces in an isodynamic 
configuration should be minimal, which is impossible, however, in view 
of the results obtained by Palumbo [11]. 

The magnetic fields that form a family of nested magnetic surfaces, 
B • n = 0, in a finite spatial region play a governing role in plasma 
confinement. Magnetic configurations for plasma confinement can be 
divided into open configurations with rippled cylindrical nested mag- 
netic surfaces and closed configurations with nested toroidal magnetic 
surfaces of complicated shape. The solenoidal nature of the magnetic 
field implies that the toroidal magnetic flux within a magnetic sur- 
face is conserved. This is why the most general equation for a family 
of nested magnetic surfaces is formulated in terms of the toroidal mag- 
netic flux: &(x,y,z) = const. The function $ is a single-valued solu- 
tion, if there is any, to the equation B- V<1? = with a known magnetic 
field having a nonzero rotational transform. Hence, at each point of 
the plasma confinement region, the vector field of unit vectors normal 
to the magnetic surfaces, n = V$/|V$| such that n • (V x n) = 0, is 
usually defined. By substituting the vector n into (2), it is possible to 
determine the mean curvature H$ of the magnetic surfaces — a quan- 
tity that plays an important role in the theory of plasma confinement 
systems. 

For completeness sake, we supplement the vectors b and n, which 
have been introduced above, with the binormal vector t = b x n up 
to an orthonormalized magnetic basis. The vector t is directed along 
the vector B x V$, which is orthogonal to B. In an equilibrium 
state in which the plasma currents j flow along the magnetic surfaces 
j • V<I> = the vector B x V<£ is solenoidal and its lines lie on the 
magnetic surfaces [9]. Applying (2) to t defines the mean curvature 
of the additional field: 

H t = t- Vlog^/B\v¥\. (6) 

This curvature is related to such familiar parameter as the geodesic 
curvature of the magnetic field lines [9] . 

2.1 Magnetic surfaces with H = const 

The energy and particle confinement in magnetic systems is commonly 
characterized by integral confinement times. For definiteness, let us 
consider the particle confinement time tn- This time is calculated 
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from the formula Tn = N/I, where N is the total number of particles 
in the system after the injection of the particle current /. In turn, the 
total number of particles is N = nV, where n is the mean particle 
density in the system and V is its volume. Since, in a steady state, 
the injection current is equal to the loss current, we have I = nvS, 
where v is the mean velocity with which the particles escape from the 
system through a boundary region of area S. As a result, the particle 
confinement time is given by the formula tn = V/vS. The better the 
confinement, the longer this time. Under the assumption that v is 
constant, we arrive at the conclusion that the boundary surfaces of 
constant mean curvature are optimum for plasma confinement. 

As shown in Section 1, there exist periodic axisymmetric undu- 
loid surfaces with H$ = const (see Fig. 2) that can be used to 
optimize plasma confinement in ambipolar open magnetic systems 
[9]. The central solenoid of such configurations with a straight mag- 
netic axis is characterized by the length L and on-axis mirror ratio 
P = B max /B m i n , or the ratio of the maximum and minimum radii 
r max/ r min = VP- The values of L and P determine the parameters of 
an ellipse that generates an optimum unduloid for a given geometry 
of the magnetic configuration. 

In Section 1, it was pointed out that tori with H$ = const do not 
exist. It is, however, for an ambipolar open magnetic confinement 
system that we consider asymptotic surfaces of constant mean curva- 
ture. Figure 3 shows a way how to close open systems by means of 
the magnetic surface of a Kadomtsev rippled confinement system [12]. 



L 




Figure 3: The Kadomtsev confinement system. 

Let us modify the Kadomtsev system as follows. On two long straight 
portions, it is possible to use unduloids with a relatively weak mean 
magnetic field, and curvilinear elements can be half-tori with a mag- 
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netic field strong enough for the radius of the tori to be small. The 
system can also be arranged to have a square shape: four straight por- 
tions with unduloids are closed by curvilinear elements in the form of 
quarter-tori. This method of modifying the Kadomtsev system can be 
repeated over and over again. The main idea is that, for a sufficiently 
large ratio L/L c , where L c is the length of the curvilinear elements, 
and a sufficiently large mirror ratio P = B c /B, where B c is the mag- 
netic field in a curvilinear element, the ratio V/S is asymptotically 
determined only by the straight portions of the system. 

2.2 Calculation of the shape of the magnetic 
Surfaces 

The main plasma confinement problem is to calculate a distortion 
of the shape of the magnetic surfaces that are in equilibrium with a 
plasma of isotropic pressure, p = p(Q). The pressure is assumed to be 
constant on the magnetic surfaces because it is equalized by a mecha- 
nism associated with the fast motion of charged plasma particles along 
magnetic field lines with irrational rotational transforms. An increase 
in pressure distorts the magnetic surfaces and, as a rule, destroys their 
nested structure, thereby deteriorating the plasma confinement rad- 
ically. This is why solving such problems is aimed at determining 
the maximum volume-averaged value of the parameter /3 eq ~ (~b?) 
that is consistent with equilibrium. This maximum beta value, and 
consequently, the plasma confinement efficiency, is very sensitive to 
the geometric shape of the magnetic surfaces. To be specific, we can 
mention innovative stellarators with a complicated three-dimensional 
geometry of helical tori — devices in which the maximum /3 eq values 
are more than an order of magnitude higher than those in tokamaks 
(see [9] for comparison) . That is why it is important to search for spe- 
cial geometries of magnetic confinement systems. In this way, it may 
be very helpful to use Eq. (2) written for the normal to a magnetic 
surface in the form of an elliptic equation for determining §(x,y,z) 
at given |V$| and H$: 



Let us introduce the function 5 = $/|V$| which is in fact the 
distance between neighboring magnetic surfaces. The function 5 so 
introduced is required to calculate the magnetic field B on a given 
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magnetic surface [13]. Instead of 5, it is possible to use the abso- 
lute value of the magnetic field B [14]. This is why the choice of the 
function |V$| is largely determined by how the absolute value of the 
magnetic field varies along its lines. Moreover, the pseudosymmetry 
condition, which is now widely used to synthesize innovative stellara- 
tors, is also constructed based on the geometry of isomagnetic surfaces 
(i.e., the surfaces of constant magnetic field strength B) [9]. 

The function Hs determines the geometry of the confinement sys- 
tem. From an expression obtained in [14], specifically, 

H S = ^^, (8) 

where K n and the normal curvatures of the lines of the magnetic 

field B and its complement B x V<£, it is clear that the pressure p 
enters Eq. (7) through the mean curvature. In fact, the (equilibrium) 
force balance equation yields the equalities 

1 „ / B 2 \ |V$| dp „, „ 
«n = l2-v( P + T J= V J + n-VlogB. (9) 

The expression for k n = n • (t • V)t takes the simplest form for ax- 
isymmetric currentless configurations [15]: 

|V$| |V$| 
k n = - ■ Vlog 1 —^- = -n- Vlog^ -n- VlogS. (10) 

Simple manipulations put Eq. (7) into the form 

V<I> dp 

(|v$| 2 /£ 2 ) = ~d$~ 

Since, in cylindrical coordinates, we have |V$| 2 /i? 2 = r 2 up to con- 
stants, we arrive at the Grad-Shafranov equation. 

Hence, after some algebraic manipulations, Eq. (7) can be used 
not only to search for an optimum geometry of the boundary magnetic 
surface but also to solve equilibrium problems. 



2.3 The integral <f s HsdS on nested magnetic 
surfaces 

Let us consider the radial variation of the integral <f s H$dS. To do 
this, we use the divergent expression presented in [2] for the Gaussian 
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curvature K of a unit vector field e: 

K= ^V- (2He + K). (12) 

where k = — e x (V x e) is the curvature vector of the lines of the 
vector field e. Setting e = n, where n is the normal to a family 
of nested magnetic surfaces and integrating the curvature from 
formula (12) over the volume dV = 4^fp between two neighboring 
magnetic surfaces $ and $ + d&, we obtain the equalities 

j> KdV = j> j^ dS = d j H sdS. (13) 

Here, we have used the fact that the magnetic flux between the surfaces 
is constant, d§ = const. The result is 

"f° H J dS _ / ^ iS . (14) 

Since § s KdS = for toroidal surfaces, and since |V<3?| ~ const, the 
integral on the right-hand side of formula (14) vanishes. 



Conclusions 

The mean curvature of the magnetic field vector is related to the vari- 
ation of the absolute value of the magnetic field along its lines. In the 
presence of magnetic surfaces and, consequently, of the orthonormal- 
ized magnetic basis (b,n, t), the mean curvature can be introduced 
for each basis vector. The mean curvature of the normal vector coin- 
cides with that of the magnetic surface. Magnetic surfaces of constant 
mean curvature, having a minimum surface area at a fixed volume, are 
optimum for plasma confinement in multimirror open systems and rip- 
pled tori with straight portions. By specifying the mean curvature of 
the magnetic surfaces and the distance to the nearest magnetic sur- 
face, it is possible to calculate the shape of the magnetic surfaces. 
All this goes to show that it may be helpful to use the notion of the 
mean curvature in the geometry of magnetic fields in plasma magnetic 
confinement systems. 
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